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We present ab initio studies of photoelectron spectra for above threshold detachment (ATD) of 
F“ anions in short, 1300 nm and 1800 nm laser pulses. We identify and assess the importance of 
electron rescattering in strong-held photodetachment of a negative ion through comparison with 
an analytic, Keldysh-type approach, demonstrating the capability of ab-initio computation in the 
challenging near-IR regime. We further assess the inhuence of the strong electron correlation on the 
photodetachment. 
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Electron rescattering is one of the fundamental pro¬ 
cesses occuring in the interaction between matter and 
intense light fields [ij. The mechanism is a critical part 
of the well known three-step or recollision model for high 
harmonic generation (HHG) or strong field double ionisa¬ 
tion. According to the model an electron is first ionised, 
then driven by a strong laser field, before recolliding with 
the parent ion, either recombining, leading to HHG ii, 
or rescattering, leading to high-energy electron emission 
ii , or non-sequential double ionisation Q. 

Electron rescattering also encodes structural informa¬ 
tion about the residual ion into the wavepacket of the 
ejected electron and can thus be exploited as an experi¬ 
mental probe of the structure of the parent ion [i[. The 
technique is especially sensitive as the current density of 
a recolliding electron wavepacket exceeds that of conven¬ 
tional electron sources by several orders of magnitude i . 
Eurthermore, the inherently subcycle and phase-locked 
nature of the recollision process gives access to electron 
dynamics on the attosecond scale, via information em¬ 
bedded in the photoelectron spectrum d, . 


One of the open questions in strong-held science con¬ 
cerns the importance of electron rescattering for nega¬ 
tive ions. Signihcant progress has been made in under¬ 
standing and controlling the equivalent process in neu¬ 


tral atoms and positive ions [10|, but above-threshold 
detachment (ATD) presents a different challenge. The 
small binding energy allows detachment at low intensi¬ 
ties. Hence to reach signihcant recollision energies, near- 
infrared (NIR) laser helds are required. In addition, the 
absence of the Goulomb potential makes it easier for the 
electron wavepacket to spread out, reducing the effect 
of rescattering ii- While evidence for rescattering 
from negative ions has been found experimentally 0, 
no verihcation has yet been provided from ab initio the¬ 
ory. A theoretical approach, based on hrst order cor¬ 
rection to the strong held approximation, was able to 
reproduce experimental results from Br“ and E“, 121 


but a more recent study, using a numerical solution of 
the time-dependent Schrodinger equation (TDSE), found 


“no qualitative evidence of rescattering” for H [13j. In 


this report we demonstrate that ab-initio theory can be 
used to investigate rescattering in the NIR regime. 

An additional complication in the description of neg¬ 
ative ions is the much larger inhuence of dielectronic- 
repulsion. Several approximate methods have been em¬ 
ployed to model photodetachment from complex negative 
ions [ISlil , but these methods are limited in their de¬ 
scription of electron repulsion. A previous study found 
that multiphoton detachment rates for E“ are affected 
substantially by the inclusion of correlation effects 0. 
The accurate description of dielectronic repulsion may 
prove especially important in rescattering calculations for 
negative ions, as the neutral core means the process will 
be mediated entirely by short-range effects rather than 
the long-range potential of a charged core. 

In this letter we address these fundamental ques¬ 
tions by applying R-matrix theory with time-dependence 
(RMT) to study strong-held dynamics of E“ in NIR laser 
pulses. The RMT method is an ab initio method for 
solving the TDSE for multielectron atomic systems in in¬ 
tense, short laser pulses. As with all R-matrix methods, 
it employs a division of space, whereby electron exchange 
effects are fully described in an inner-region close to the 
nucleus, while far from the nucleus, a single, ejected elec¬ 
tron moves in the long-range multipole potential of the 
core. Although several other time-dependent R-matrix 
methods have emerged in recent years [iM3, RMT 
exhibits orders-of-magnitude improvements in efhciency, 
primarily because it employs hnite-diherence (ED) tech¬ 
niques to model the one-electron wave function in the 
outer region. RMT merges the outer region ED model 
with a B-Spline-based, R-matrix basis set for the mul¬ 
tielectron inner region, allowing efficient calculations ac¬ 
counting for atomic structure and correlation effects 0 . 

The R-matrix basis for E“ consists of the neutral E 
ground state coupled with an additional electron. We 
employ two different models for the E atom. These two 
models are compared to assess the influence of electron 
correlation. The basic model includes only one configura¬ 
tion, ls^2s^2p^, with the Is, 2s and 2p orbitals given by 
the Hartree-Eock orbitals for the E ground state [15|, |20j . 
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FIG. 1: (Color online) Electron momentum distributions for ATD of F~ in the — kz plane for a (a) 1300 nm (b/c) 1800 nm 
laser pulse of intensity 1.3 x 10^® W/cm^. (a) and (b) are calculated using the RMT approach, while (c) is obtained using the 
KTA method [ 13 ]. 


The second model includes additional 3s, 3p, and Zd pseu¬ 
doorbitals, [2l|. This allows us to generate an accu¬ 
rate wave function expansion for the ls^2s^2p® ^P° state 
from a configuration interaction calculation including the 
ls^2s^2p®, ls^2s2p®3s, ls^2s^2p"^3p, ls^2s^2p^3p^ and 
\s^2s^2p^Zd?‘ configurations. Model 2 gives a binding en¬ 
ergy of 3.421 eV for the ^5'® F ground-state, which is close 
to the experimental value of 3.399 eV [ 23 , whereas in 
model 1 we shift the ground-state artificially to 3.421 eV. 

We employ 1300 nm and 1800 nm wavelength fields 
at a peak intensity of 1.3 x 10^^ Wcm“^. The profile 
comprises two cycles sin^ ramp-on followed by two cycles 
sin^ ramp-off. In such high intensity, long-wavelength 
fields, the ejected electron wavepacket can travel far from 
the nucleus, and hence, an accurate description of the 
wave function is required over an extended region of 
space. To facilitate this we include angular momenta 
up to L = 240 and propagate the wave function out to 
a radius of 4265 ag. The time step in the wave function 
propagation is 0.24 as. 

The laser parameters are chosen to facilitate a com¬ 
parison with results obtained using a Keldysh-Type Ap¬ 


proach (KTA) calculation [IJ, [23|, [2^ . In this model, the 


effect of the atomic potential on a detached electron is ne¬ 
glected. The description of F“ is based on Hartree-Fock 
orbitals for F“ in which the long range tail is fitted to the 
correct binding energy. The laser field is infinitely long, 
but with a periodic envelope which describes a series of 
short pulses. This allows the analytical solution of the 
so-called saddle-point equation, yielding the electron tra¬ 
jectories in the field. As a consequence of the long laser 
‘pulse’, the field is composed of three distinct photon en¬ 
ergies, which is at variance with our time-dependent R- 


matrix calculations wherein a spread of photon energies 
follows from the isolated short pulse used. 

Figure [T] shows the 2D electron momentum spectra for 
ATD from F“ for 1300 nm and 1800 nm fields alongside 
the KTA spectrum for 1800 nm 141. The z-axis is di¬ 
rected along the polarization axis of the laser. The figure 
shows results from the basic atomic structure, model 1 of 
F“ described above, which allows a closer comparison to 
the KTA results. The F atom has three different possible 
final states, corresponding to m = 0, ±1 where m is the 
magnetic quantum number. The total momentum distri¬ 
bution is obtained by incoherently summing the contri¬ 
butions from the m = 0 and the m = ±1 states. 

The ATD momentum distribution in Fig. [I] is exten¬ 
sively detailed. Both the RMT and KTA spectra show 
the typical structure of rings centered at zero momen¬ 
tum, with each ring corresponding to the absorption of 
N photons. Interference patterns arise because multiple 
electron trajectories contribute to each final momentum 
state 13|. Lines connecting the interference minima are 


curves in momentum space which satisfy a destructive 
interference condition. Curves intersecting the P^-axis 
correspond to interference between the well known long- 
and short-trajectories. Curves intersecting the P^j-axis 
lead to minima/maxima in the even/odd ATD rings for 
perpendicular emission. The two different approaches 
therefore demonstrate the same basic physics of photode¬ 
tachment. However, the features appear sharper in the 
KTA spectra due to the well defined photon energy, w, 
while for RMT the results are blurred by uncertainty in 
w. Furthermore, the RMT spectrum appears narrower in 
the kx direction. 

Figure [T] shows additional ATD rings for electron mo- 
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menta greater than the maximum drift momentum, Pz = 
Eo/oj (where Sq is the maximum strength of the electric 
field) for both 1300 nm and 1800 nm with the RMT cal¬ 
culations. These extra rings are not present in the KTA 
model results. In the positive z direction, the rings are 
centred on Pz sa 0.55 a.u., and Pz ^ 0.75 a.u. for 1300 nm 
and 1800 nm respectively: these values corresponding to 
the maximum drift momentum for each wavelength. 

These additional rings are a clear signature of electron¬ 
rescattering. In the strong field electrons are detached 
near the maximum of the held each half cycle as per the 
three-step model. The electrons are thus ‘born’ into the 
held with zero energy, leading to ATD rings centered on 
zero momentum as in Fig. [TJ Recollision then occurs near 
the zero-held intersections, allowing the electron to gain 
extra drift momentum through its interaction with the 
held potential. Rescattered electrons are effectively ‘re¬ 
born’ into the held at the held-zero, gaining a maximum 
drift momentum of Pz = EojoJ- This can be conhrmed 
with classical trajectory calculations [1^. The appear¬ 
ance of these rings is proof of the recollision mechanism 
in negative ions. Angle resolved spectra showing these 
rings centered on the drift momentum go beyond the in¬ 
tegrated yields shown in Fig. [2] and correspondingly in 
[l2l |. which show only evidence of high energy electrons 
which are ascribed to the recollision mechanism. 

The results display a clear asymmetry in the Pz direc¬ 
tion, with the recollision rings in the positive z direction 
extending to larger momenta than in the negative z di¬ 
rection. This is due to the short pulse profile we employ. 
The highest energy electrons are accelerated by the sin¬ 
gle, peak-intensity cycle of the laser pulse and hence are 
only emitted in the positive direction. Moreover, in the 
negative z direction, there is a pronounced interference 
structure in the recollision rings, arising from two inter¬ 
fering electron trajectories in the ramp-on and ramp-off 
of the field. Finally, the centre of these recollision rings is 
shifted towards the origin for the negative z direction, be¬ 
cause the electron trajectories in the ramp-on and ramp- 
off of the field are driven by lower intensity laser cycles. 
Thus the centre of the recollision rings in the negative 
z direction is given by 0.85 x SqIuj. We note that the 
non-recollision part of the momentum distribution is me¬ 
diated by the vector potential of the laser field, which is 
antisymmetric about zero potential. Thus this part of 
the distribution is symmetric with respect to the origin, 
while the recollision, mediated by the asymmetric elec¬ 
tric held, leads to an asymmetric spectrum. We would 
expect the opposite to be true, were the carrier envelope 
phase of the laser pulse adjusted accordingly. 

A second outstanding question surrounding strong- 
held dynamics of negative ions concerns the inhuence of 
electron correlation. Since the outer electron in F“ is 
loosely bound, and the main rescattering process occurs 
very close to the residual atom, electron correlation may 
affect both the initial ionization and the rescattering. A 


TABLE I: The level of ionisation of F“ in 1300 nm and 
1800 nm fields calculated for different peak intensities. This 
ionised population is presented for the two models used in 
RMT, compared to the total detachment probability from 0 


Intensity 


Detachment probability 


W/cm^ 

RMT(Model 1) RMT(Model 2) 

Ref. [14] 

1300 nm 

7.7 X 10^^ 

0.018 

0.011 

0.036 

1.1 X 10^® 

0.045 

0.031 

0.090 

1.3 X 10^® 

0.065 

0.044 

0.139 

1800 nm 

7.7 X 10^^ 

0.020 

0.013 

0.023 

1.1 X lO^^^ 

0.055 

0.034 

0.066 

1.3 X 10^® 

0.080 

0.052 

0.106 


recent study suggests that while for neutral atoms there 
are notable differences between KTA model results and 
those from the numerical solution of the TDSE, these 
differences are negligible for negative ions [l3|. The 
ejected-electron momenta spectra obtained for model 1 
and model 2 show little qualitative difference, but no¬ 
table, quantitative differences can be observed when elec¬ 
tron correlation is included. 

Table U gives the population in the outer region— a 
measure of detachment probability in our calculations— 
compared with the detachment probability from the KTA 
model for 1300 nm and 1800 nm. The yield in our model 

1 calculations is reduced by a factor of 2 from the KTA 
model result for 1300 nm, whereas the difference is only 
about 20% for 1800 nm. When using the more sophis¬ 
ticated model 2 for F“, we find a further reduction in 
the detachment probability by one third. The reason for 
these differences is not obvious, although the effective po¬ 
tentials differ. The binding energy in RMT model 1 is 
shifted artificially to the experimental value, while model 

2 gives the correct binding energy directly. This implies 
that the short-range potential, and therefore the wave 
function, is described more accurately in model 2. This 
increase in accuracy may lead to detachment yields about 
a third smaller than model 1, and a factor of 2-3 smaller 
than the KTA. Differences in short-range potential can 
thus lead to significant differences in detachment yield. 

The differences between model 1 and model 2 are most 
easily identified in the photelectron energy spectra, which 
are presented in Figs. [2] and |3] alongside results from the 
KTA model. In all calculations, model 1 shows a larger 
electron yield across the full energy range, reflecting the 
lack of change in the overall structure of the spectrum. 
The figures show that the main effect of electron interac¬ 
tion is an overall reduction in detachment probability. 

Figure [2] demonstrates the capability of the RMT ap¬ 
proach to describe electron rescattering. For energies up 
to about 10 eV, the RMT and KTA calculations show a 
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FIG. 2: (Color online) Photoelectron energy spectra for ATD 
from F“ by a 1.3 X 10^® W/cm^ pulse at 1300 nm and 1800 nm 
from the present approach, model 1 (solid, red line), model 
2 (dashed, blue line) and the KTA model approach (dotted, 
black line). The presence of ‘knee’-like structures in the RMT 
results is a characteristic signature of rescattering. 


very similar behaviour for electron emission, consisting of 
a small initial plateau followed by rapid exponential de¬ 
cay. Differences emerge for energies greater than about 
10 eV, with the RMT calculations showing a plateau at a 
magnitude of 10“^-10“® eV“^ extending to an energy 
of about 20 eV for 1300 nm and 35 eV for 1800 nm, cor¬ 
responding to around lOC/p in each case. This plateau is 
absent in the KTA calculations. Since the KTA approach 
does not account for rescattering we can identify this as 
the source of the plateau. The yields obtained in model 
1 and model 2 are very similar, with those obtained in 
model 2 lying about one third lower than for model 1. 

Figure [3] shows the low energy region of the photo- 
electron spectra, which display peaks characteristic of 
the multiphoton mechanism: each peak corresponding 
to the absorption of an integer number of photons and 
with energy, E = Nui — Up — Ip for some integer N, and 
the ionization energy. Ip. In all cases, the magnitude of 
the photoelectron peaks is lower for the present calcula¬ 
tions than for the KTA approach and lower for model 2 
than for model 1, consistent with the lower detachment 
probability discussed above. The qualitative agreement 
with the KTA model is excellent for 1800 nm, while, at 
1300 nm, the positions of the ATD peaks differ. The 
origin of this difference is difficult to assess unambigu¬ 
ously. However, two factors may play a role. Firstly, the 
RMT calculations account for AC-Stark shifts beyond the 
ponderomotive energy shift, whereas the KTA model in¬ 
cludes the ponderomotive shift only. Secondly, the RMT 
method uses a short laser pulse, while the KTA uses an 
infinitely long field with a periodic envelope. This leads 
to differences in the photon energy and consequently the 
ponderomotive shift, making it difficult to determine the 
correct position of the photoelectron peaks. It will there¬ 
fore be interesting to compare these results with exper¬ 
imental spectra, although, as far as we are aware, only 


FIG. 3: (Color online) As Fig. [2]but shown for the low energy 
region on a linear scale. 


on-axis emission has been measured for the present laser 
parameters [2^ . 

In conclusion, we have demonstrated the capabilitiy of 
ab initio theory to study ATD of F“ anions in computa¬ 
tionally challenging NIR laser pulses with a full descrip¬ 
tion of multielectron effects. Through comparison with a 
KTA model we have identified the importance of the rec¬ 
ollision mechanism in the electron momentum distribu¬ 
tions. High energy rings in the angle resolved photoelec¬ 
tron spectra are the first theoretical verification of rescat¬ 
tering from negative ions. The integrated photoelectron 
energy spectra shows further clear evidence of rescatter¬ 
ing, with a knee-like, plateau structure extending to en¬ 
ergies of 10C/p. Although the yield in these rescattering 
channels is small- on the order of .02 % and .07 % of 
the total yield for 1300nm and 1800nm respectively- it 
is clear evidence of rescattering. 


The description of multielectron atoms and ions in 
short, intense NIR pulses is a significant computational 
challenge, but holds substantial promise for elucidat¬ 
ing new physics, mediated by multielectron interaction. 
The RMT method has already been applied to other 
strong held phenomena, including high-harmonic genera¬ 


tion 23 , but the accurate and efficient determination of 
the multielectron wave function over such a large region 
of space is a first for calculations of this type. Indeed, 
the present calculation of ATD in negative ions provides 
an extremely sensitive test of the wave function accu¬ 
racy. It will be interesting now to compare theoretical 
single-atom results of the RMT approach with the find¬ 
ings of state-of-the-art experimental techniques in strong- 
held physics and model approaches 0. 
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The data used in this paper may be accessed via 
pure.qub.ac.uk/portal/en/datasets/search.html 
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